Additive-multiplicative stochastic models of financial mean-reverting processes 
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We investigate a generalized stochastic model with the property known as mean reversion, that 
is, the tendency to relax towards a historical reference level. Besides this property, the dynamics 
is driven by multiplicative and additive Wiener processes. While the former is modulated by the 
internal behavior of the system, the latter is purely exogenous. We focus on the stochastic dynamics 
of volatilities, but our model may also be suitable for other financial random variables exhibiting the 
mean reversion property. The generalized model contains, as particular cases, many early approaches 
in the literature of volatilities or, more generally, of mean-reverting financial processes. We analyze 
the long-time probability density function associated to the model defined through a Ito-Langevin 
equation. We obtain a rich spectrum of shapes for the probability function according to the model 
parameters. We show that additive-multiplicative processes provide realistic models to describe 
empirical distributions, for the whole range of data. 

PACS numbers: 89.65.Gh, 02.50.Ey, 05.10.Gg 



I. INTRODUCTION 

Accurate statistical description of the stochastic dy- 
namics of stock prices is fundamental to investment, op- 
tion pricing and risk management. In particular, a rel- 
evant quantity is the volatility of price time series 
that quantifies the propensity of the market to fluctuate. 
Since volatility represents a measure of the risk associ- 
ated to the fluctuating dynamics of prices, it is crucial to 
develop suitable models to predict its complex intermit- 
tent behavior. There is empirical evidence that it fluc- 
tuates following a stochastic dynamics subjacent to that 
of prices, whose dynamics, in turn, depends on the time 
evolving volatility. Many approaches are based on that 
assumption 2j, although others propose the existence of 
a reciprocal feedback between both processes Q. 

Our approach builds on the development of a simple 
Langevin equation to characterize the stochastic process 
of volatility. The equation provides an unifying descrip- 
tion that generalizes widely discussed models in the liter- 
ature. We analyze the shape of the long-time probability 
density function (PDF) associated to the stochastic dif- 
ferential equation that characterizes each particular case 
of the generalized model. Most previous results focus on 
the tails of the PDFs. In fact, for stochastic variables, 
such as volatilities presenting fat-tailed PDFs0,O|, it is 
specially important to reproduce extreme events in a real- 
istic model. Now we go a step further and aim to predict 
the PDFs in the whole range of events. 

One of the main features observed in the dynamics of 
some financial variables, such as volatilities, stock vol- 
umes or interest rates, is their tendency to permanently 
relax, towards a reference level 9, a property known as 
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mean reversion. Another feature is the multiplicative 
market processing of random news, whose strength be- 
comes modulated by a function of the stochastic variable 
itself. These two properties are modeled by means of a 
nonlinear mean-reverting force and nonlinear multiplica- 
tive noise. They are discussed in detail in Sect. [H] 

In Sect. IIIII we discuss the shapes of the PDFs that 
such family of models yields. Despite being of a gen- 
eral form, they give rise to PDFs that, decay exponen- 
tially fast, either above the mode, below it, or both, in 
disagreement with empirical observations. For instance, 
log-normal behavior has been reported for volatility com- 
puted from global data of S&P5OO0, at intermediate val- 
ues. However, at high values, a power-law behavior, with 
exponent outside the stable Levy range, was observed. 
The same analysis performed for individual companies Q 
yields also power-law tails. But in that case, the results 
show a variation slower than log-normal below the mode, 
suggesting a power-law also in the limit of small values. 
The volatility of capitalized stocks traded in US equity 
markets exhibits similar features^. Other variables with 
mean-reversion, such as volume of transactions (number 
of trades) present akin distributions. Power-law tails out 
of the Levy range have been reported for the PDFs of 
normalized NYSE stock volumes J\ . More recently, stud- 
ies of normalized volumes, performed over high resolu- 
tion data (1-3 minutes) of NYSE and NASDAQ d (see 
also |j), display PDFs with power-law behavior both at 
large and small values. We will show that the class of 
multiplicative processes considered in Sect. IIIII although 
general enough, is not able to reproduce, for any value of 
its parameters, these empirical PDFs in the whole range. 

In a realistic model, we must deal with various sources 
of fluctuations acting upon the collective variable. Then, 
we propose to include a component that is lacking to suit- 
ably model many real processes, that is the presence of 
fluctuations that act additively, besides the multiplica- 
tive noise already taken into account. The latter origi- 
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nates from the internal correlated behavior of the mar- 
ket, representing a sort of endogenous feed-back effect, 
while additive noise concerns fluctuations of purely exter- 
nal origin or random speculative trading. Then, in Sect. 
IIVI we present a further generalization that consists in 
incorporating an independent additive source of noise. 
Depending on the parameters of the process, the addi- 
tive or multiplicative contributions will play the domi- 
nant role. This gives rise to a rich spectrum of PDF 
shapes, in particular, a subclass with two-fold power-law 
behavior, both above and below the mode, providing a 
general realistic framework for describing the shape of 
empirical distributions. A comparison with experimen- 
tal results is presented in Sect. E] Finally, Sect. IVfl 
contains the main conclusions and general remarks. 



II. MEAN REVERSION AND 
MULTIPLICATIVE FLUCTUATIONS 

The reversion to the mean is one of the basic ingredi- 
ents to describe the dynamics of several stochastic vari- 
ables of interest in economy, ft is fundamental since it 
concerns the behavior around a central value 6 and re- 
flects the global market response to deviations from a 
consensus or equilibrium level. It depends on monetary 
unit, market size, degree of risk aversion, etc., hence, it is 
characteristic of each market. The aversion to deviations 
from the mean needs not be linear, specially when large 
deviations are involved. Similarly, a nonlinear mecha- 
nism due to the cooperative behavior of traders, rules 
the way the market modulates the amplitude of fluctua- 
tions (mainly external) giving rise to innovations. 

We consider the general class of stochastic differential 
equations given by 

dx = -j[x - 9]x r ~ 1 dt + px s dw , (1) 

where, 0, 7, > 0, r, s G R, and w is a Wiener process, 
such that (dw) = and ({dw) 2 ) = 2dt. The defini- 
tion of the stochastic process is completed by the ltd 
prescription. This class generalizes well-known models 
employed to describe the dynamics of mean-reverting fi- 
nancial variables 10]. In particular, some traditional pro- 
cesses for modeling volatilities or, mainly, squared volatil- 
ities are the Hull- White (r = 1, s = l)[ll| and the Heston 
(r = l,s = 1/2)0] models, the latter also known ei- 
ther as Cox-Ingersoll-Ross fl3j or Feller process 0j. The 
arithmetic (r = 1,8 = 0) and geometric (r = 2, s = 1) 
Ornstein-Ulhenbeck processes are particular cases too. 
Moreover, several other models employed in the litera- 
ture of volatilities are related to this class |rj fill. 

Different values of r in Eq. represent different pos- 
sible relaxation mechanisms of amplitude 7, determined, 
amongst other factors, by constraints, flux of informa- 
tion, stock liquidity and risk aversion, which are partic- 
ular of a given market. Notice that the restoring force 
in Eq. corresponds to a confining potential, with 



minimum at 9, for all r £ R. The larger r, the more at- 
tractive the potential for large x, but the less attractive 
for vanishing x. Similarly, different values of s specify 
the market informational connectivity, which conditions 
the degree of coherent multiplicative behavior. Models in 
the literature typically set s > 0, meaning that the effec- 
tive amplitude of fluctuations increases with x. Negative 
s makes multiplicative fluctuations grow with decreasing 
x, thus it mainly reflects a cooperative reaction to qui- 
escence. Although it does not seem to reflect a realistic 
steady state of the market, it may occur as a transient, 
driven by speculative trading. 

The two mechanisms are complementary. If r < the 
restoring force decreases for increasing x above the ref- 
erence level, in particular, for r < — 1, the force tends 
to zero in the limit of large x. Thus, decreasing r repre- 
sents markets that, become less able to recover the refer- 
ence level by means of the deterministic tendency alone. 
However, a strong multiplicative response to large fluc- 
tuations (positive s) could still compensate that inability 
and restore the market historical level. Concerning the 
response to small values, the restoring force diverges at 
the origin if r < 1, while for r > 1, it vanishes at x = 0, 
meaning that this point becomes an unstable equilibrium 
state. This corresponds to a market indifferent to low lev- 
els of trading activity. Again, this effect can be balanced 
by the multiplicative contribution (with a small value of 
parameter s). 

In early works, only very particular values of (r, s) 
have been considered. However, this may be sometimes 
owed more to reasons of mathematical solvability, than 
to econophysical ones. Following the above discussion, 
(r, s) may be non-universal, depending on the particular 
nature of a market or its agents. Therefore, we will not 
discard any possibility a priori. 

III. GENERALIZED MULTIPLICATIVE 
PROCESS WITH MEAN REVERSION 

We consider the simple class of stochastic multiplica- 
tive differential equations given by Eq. 0, that gener- 
alizes many processes usually found in the literature of 
volatilities. We investigate, in this Section, the long-time 
PDFs that this class of processes yields. The Fokker- 
Planck equation associated to Eq. QJ , following standard 
methods [lq. is 

dtp = ld x {[x - 9]x r - l p) + p 2 d 2 xx [x 2s p] . (2) 

Its long-term solution is relevant in connection to the 
assumption that the process can be treated as quasi- 
stationary. In that case the PDF obtained from an ac- 
tual data series will coincide with the stationary solution. 
Considering reflecting boundary conditions at x = and 
x — > oo[lq|. the steady state solution of Eq. reads: 

x zs \ lp + 1 p 1 J 
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with p = r — 2s ^ 0,-1, where p is a normalization 
constant and 7 M = 7//i 2 an effective restoring amplitude, 
such that 7 (a parameter associated to order) becomes 
reduced by the amplitude of multiplicative noise (associ- 
ated to disorder). 
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FIG. 1: Diagram of the asymptotic behavior of the PDF given 
by Eq. J3J, in (s,r)-plane. Unshadowed regions and dotted 
borders identify regions excluded by the normalization condi- 
tion. At the positive r-axis, the PDF is finite at the origin. 
Tilted lines denote the marginal cases r = 2s (p = 0), with 
pure exponential tail and power-law growth at the origin, and 
r = 2s — 1 (p = — 1), with power-law tail and exponential of 
1/x growth at the origin (the threshold points of these lines 
have coordinates ([1 + r y ll 6]/2, 0) and ([1 — 7 M ]/2, 0), respec- 
tively). Parameter a > 0, in the exponential formulas, as 
well as the power-law exponents, depend on model parame- 
ters. Symbols correspond to the special processes: Hull- White 
(HW), Heston (H), Ornstein-Ulhenbeck (OU) and geometric 
OU (GOU). 

The class of processes described by Eq. @ thus gener- 
ically yields asymptotic exponential-like behaviors for 
small or /and large values. As soon as p + 1 > 0, a 
stretched exponential decay is obtained for large enough 
x, such that the argument of the exponential is propor- 
tional to —x p+1 . If p < 0, a stretched exponential of the 
inverse argument (— I/a;' p ') is obtained for vanishing x. 
Therefore, for p £ (—1,0), the PDF presents dominant 
exponential-like behavior both for low and large values, 
without any restriction on the value of s. Outside that 
interval, the power law x~ 2s in @ asymptotically dom- 
inates, for either small (if p > 0) or large (if p < — 1) 
argument. Then, normalization in [0, oo) restricts the 
possible values of s according to: s < 1/2 (if p > 0), 
x > 1/2 (if>< -1). 

In the marginal cases, Eq. Q explicitly is: 

I: For p = i — 2s = — 1, 



p(x) 



Po 



exp(-7 M 6>/a;) 



(4) 



with 2s > 1 — 7 M for nomalizability. 
II: For p = r - 2s = 0, 

p(x) = j o x 7 " e - 2s exp(-7 M x) 



(5) 



with 2s < 7^0 + 1 for normalization, but 2s < j^O to 
avoid the divergence at the origin. 

Fig. ^ displays the possible PDF asymptotic shapes 
in (s,r) space. Notice that the s = axis gives the 
solution for mean-reverting models with purely additive 
fluctuations. Let us analyze some special cases. In the 
trivial case (s, r) = (0, 1), corresponding to the Ornstein- 
Ulhenbeck process [Ty| 



d.r 



-^[x — 6]dt + pdz 



(6) 



the noisy contribution becomes additive and the station- 
ary PDF is Gaussian (truncated at x = 0). 

Although we are dealing with > 0, it is worth of 
mention the case (0, s, r) = (0, 1, 1), corresponding to the 
geometric Brownian process, that leads to the log-normal 
distribution. 

For type I (r = 2s — 1), notice that the PDF decays 
as a power law, for large x, and goes to zero faster than 
power law, for vanishing x. The power-law exponent is 
controlled by s and 7 M , that is, all the model parameters, 
except 9 are involved. In the particular case r = s = 1, 
one recovers the Hull- White process [Tlj 



da; = — 7 [a; — 9]dt + pxdz . 



(7) 



In case II (r = 2s), observe that the PDF has opposite 
behavior: it increases at the origin as a power law and 
decays exponentially for large x. All the model parame- 
ters, including 9 are buried in the power-law exponent. In 
particular, if r = 2s = 1, one gets the Heston model |T^ 



d:r 



-7F- 



9]dt + p^/xdz 



(8) 



If r = 2s = 2, the geometric Ornstein-Uhlenbeck process 
is obtained 



da; = 



-7F 



cdt + pxdz 



(9) 



Diverse other models proposed in the literature can 
also be thought as particular instances of our generalized 
model. For example, the one proposed by Micciche et al. 
is in correspondence with the Hull- White model JJJ), 
with x representing volatility v, whereas in the latter 
x = v 2 . Also a family of multiplicative models, stud- 
ied before in the context of a wide spectrum of phys- 
ical processes [T^. belongs to the class here considered, 
through the transformation x — > x@ . 

Summarizing, from Eqs. (J3J)-©, in general, the asymp- 
totic behaviors below and above the mode are tied, such 
that, in a log-log scale, if one flattens the other changes 
rapidly. This explains why models of this class fail to 
describe empirical volatilities in the whole range of ob- 
served data, even under the transformation v 2 i— > v. 
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IV. GENERALIZED MODEL WITH 
ADDITIVE-MULTIPLICATIVE STRUCTURE 

We analyze in this section, processes that take into 
account the presence of some additional source of 
noise. Previous works[H [H HI show that additive- 
multiplicative stochastic processes constitute an ubiqui- 
tous mechanism leading to fat-tailed distributions and 
correlated sequences. This extra noise represents a quite 
realistic feature, since, besides noise modulated by the 
market, other fluctuations may act directly, additively. 
From the stream of news, represented by a noisy signal, 
some are amplified or reduced by cooperative actions, 
others incorporated unaltered. Related ideas has been 
discussed in Ref. j2l|. Also, a model of financial markets 
that leads to additive-linear-multiplicative processes has 
been recently proposed[22j, where the noises are identi- 
fied with the fluctuating environment and fluctuating in- 
teraction network, respectively. In general, the two white 
noises are considered uncorrelated. However, they may 
even correspond to identical time-series as soon as they 
are shifted with a time lag greater than the correlation 
time. In such case, the endogenous noise is expected 
to act with a delay due to its very nature of feedback 
process, whereas, the additive noise is incorporated im- 
mediately, free of signal processing. 

By including purely exogenous fluctuations, in the pro- 
cess defined by Eq. QJ, we obtain the following Ito- 
Langevin equation (ILE) 



dx = 



-7F- 



l dt + p,x s dw + adz 



(10) 



where w, z are two independent standard Wiener pro- 
cesses, defined as above, and p, a their respective ampli- 
tudes. The corresponding Fokker-Planck equation reads 



dtp = ld x {[x - 6]x r ' l P ) + d xx ([p 2 x 2s + a» . (11) 

Its steady state solution with reflecting boundary condi- 
tions is 



p{x) 



Po 



l + A 2 a; 2seXp[ la j 1 + AV 



y r - l {y-e) 



dy], (12) 



with p a normalization constant, j a = 7/a 2 , A 2 = 
(p/a) 2 = ja/jfi- In most cases the integral can be writ- 
ten in terms of hypergeometric functions 2-F1 [2j| , through 



x y"- 1 
1 + AV< 



■dy= 4 2 Fi(c,1,c+1,-A 2 .t 2s ) (13) 



with c = /3/(2s) 7^ —1, —2, . . ., whereas, in the marginal 
case (3 = 0, we will use 



x y- 1 



1 + X 2 y 



2, ,2s 



dy ee lnx - ln(l + A 2 x 2s )/(2s). (14) 



By means of these definitions and their asymptotic 
formulas p3l |24| . we obtain the possible PDF shapes, 



in (s, r)-space, as schematized in Fig. |2 The marginal 
cases r = and r = — 1 will be considered latter. In 
general, sufficiently large positive s is required in order 
to yield power-law tails, otherwise, stretched exponential 
tails prevail, as for the processes considered in Sect. II I II 
The additive noise does not add new domains with power- 
law tails, although regions with stretched exponential law 
are excluded or included by the normalization condition. 
For vanishing x, the main difference with purely mul- 
tiplicative processes is that, for positive both s and r, 
the PDF is truncated at the origin. Notice that, as the 
PDF is finite at the origin, then, if x is identified with 
the squared volatility (x ee v 2 ), the PDF for v increases 
linearly at the origin. 

Let us analyze, in more detail, the marginal cases r = 
and r = — 1 that can yield power-laws in both asymptotic 
limits. From Eqs. ijT^ - lfHfl) . we obtain 

A: For r = 0, the PDF has the form 



p{x) = p 



Q(x) 



[1 + X2 x 1syf a 6/{2s) + l 



(15) 



where 0(x) = exp[— 7^x2-^1(^,1,^+1, — A 2 ir 2s )] is a 
smooth function of x, such that 0(0) is finite, hence it 
does not spoil the power-law growth at the origin. For 
large x, it may present different asymptotic behaviors 
depending on the value of s: 

A.l: If s < 0, &(x) decays as pure exponential of x. 
Therefore, the asymptotic decay is finally dominated by 
this exponential factor. 

A. 2: If < s < 1/2, Q(x) behaves asymptotically as 
a stretched exponential with argument x 1_2s . That is, 
the tail, although a power law for moderate x, becomes 
asymptotically dominated by a stretched exponential de- 
cay 

A. 3: If s > 1/2, Q(x) tends to a positive value, there- 
fore, in this instance, the tail remains power-law. 

There, by switching s, one tunes the tail type, being 
a power-law for s > 1/2. In the threshold case s = 1/2, 
we have 2-Fi(l, 1, 2, —z) = ln(l + z)/z, then we get the 
explicit expression 



p A (x) = Po 



,7c 6 



[1 + \2x}^ 6 +^ + 1 



(16) 



Thus, the case r = 0, s > 1/2 allows one to model 
empirical PDFs with twofold power-law behavior. 

B: In the case r = — 1, the normalization condition 
requires: s > 1/2, or also, if 7 a > 1, s < is allowed. 
The PDF has the form 



P(x) = Pol 



"7. 



9(x) 



r 2s]- la /(2s) + l ' 



(17) 



where, 6(x) ee exp[— j a 9 2 Fi(-^, 1,-^+1, -\ 2 x 2s )/x] 
tends to a finite value for large x, therefore, the tail is a 
power-law. The asymptotic behavior of Q(x) for small x, 
depends on s. 

B.l: For s > 1/2, it behaves as an exponential of — 1 / 'x , 
that dominates the low x behavior of the PDF. 



5 



B.2: For — 1/2 < s < 0, <d(x) behaves as an exponen- 
tial of —l/x 1+2s , that dominates the asymptotic behav- 
ior. 

B.3: However, 9(x) takes asymptotically a finite 
value, if s < —1/2; hence, the complete expression in- 
creases at the origin as a power-law. 

At the threshold value s = —1/2, by employing again 
the explicit expression for 2-Fi(l, 1, 2, —z), one obtains 

PB{X) = ^ [l + s/A^+^+i ' if 7q > 1 ■ (18) 

Thus, the case r = —1, s < —1/2 also provides twofold 
power-law distributions. 



the extra additive noise, new types of asymptotic behav- 
iors emerge. Specially interesting solutions arise in the 
marginal cases r = 0, — 1 where two- fold power-law PDFs 
are found. 

Moreover, additive-multiplicative processes lead to 
higher richness of crossover behaviors, with respect to 
purely multiplicative processes. Therefore, the appear- 
ance of new PDF shapes exceeds the one resulting from 
the mere analysis of the asymptotic regimes. This is spe- 
cially important because depending on the values of the 
parameters, the true asymptotic regime might fall out- 
side the observable range. 

V. COMPARISON WITH EMPIRICAL 
DISTRIBUTIONS 



exp(-ax p+1 )/(x- 
finite i'{xsO)a 




exp(-ax' +1 ) (x 
power law (x 



power law {x->°°) 
finite (x->0) 



power law (x->oo) 
exp(-a/jc IH ) (x->0) 
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FIG. 2: Diagram of the asymptotic behavior of the PDF de- 
fined by Eq. 11211 . in (s,r)-plane. Unshadowed regions and 
dotted borders are regions excluded by the normalization re- 
quirement. At both positive semi-axes, the growth at the 
origin is power law. On dark gray lines, tails are power law, 
with the tilted line corresponding to r = 2s — 1 (p = —1). 
Dashed lines correspond to pure exponential tails, with the 
tilted line corresponding to r — 2s (p — 0). In the formu- 
las, a > 0, as well as the power-law exponents, generically 
depend on model parameters, moreover p = r — 2s. Symbols 
correspond to the special processes A [Eq. I|16fl] and B [Eq. 

CHI]. 

In general, the class of asymptotic behavior is ruled by 
(s, r) that determine the form of market laws. This holds, 
of course, as soon as the remaining parameters assume 
moderate values. For instance, the factor A 2 = 70/7^ 
accompanies x 2s in the formula for p(x) [Eqs. H12fl - I|14|) ]. 
then, extreme values of A will change the asymptotic 
regime. In fact, in the limit a — (negligible additive 
noise, corresponding to A — > oo), different laws arise, as 
we have seen in the precedent Section. 

Summarizing, we have shown the whole picture of 
asymptotic behaviors that a general class of additive- 
multiplicative processes produce. As a consequence of 



Let us consider, as paradigm of the PDFs with two- 
fold power-law behavior, Eqs. (|16|) and l|18|) . that have 
a simple exact expression. Actually they have the same 
functional form, via redefinition of parameters (9, j a , 7^). 
This expression has been recently proposed in the liter- 
ature as an ansatz for fitting the distribution of high- 
frequency stock- volumes 8] , under the form 



P(x) = po- 



(X/X Y 



[1 + (q - l)x/x ] 



(19) 



where, in that specific application, x is identified with 
normalized stock volume. Therefore, identification of the 
process for real volumes with one of the models above, 
may allow an econophysical interpretation of the fitting 
parameters. Table I presents the correspondence between 
the parameters of Eq. (|19|l and those of processes A and 
B, given by Eqs. (|Ttj)l and (|T%1) . respectively. 



Eq. Ipjl 


A 

Eq. Ipjl 


B 

Eq. lO 


V(?-l) 


l + 7<^ + 7m 


1 + 7c + 7m# 


V 


7c 


1 + 7m# 




(3 - 1)7m/7q 


(<7 ~ 1)7<*/7m 



TABLE I: Correspondence amongst model parameters. 



Recall that j a = 7/a 2 and 7^ = 7/^i 2 (A 2 = j a /~fn), 
thus, the power-law exponent for small values of x, given 
by v (see Table), increases with 7 and 9, and is reduced 
by either one of the two noise amplitudes: the additive 
noise in process A and the multiplicative one in process 
B. The power-law decay (with exponent l/(q — 1) — v) 
for large values of x is ruled by either one of the effec- 
tive coefficients 7^ (in A) or -f a (in B) (see Table). That 
is, the tail is fatter, the larger the corresponding noise 
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amplitude. While in process A the multiplicative noise 
affects the tail, in model B it is affected by the additive 
noise, oppositely to what happens for small values. This 
is related to the sign of s, indicating higher multiplica- 
tive feedback for either increasing (A) or decreasing (B) 
values of x. 

Besides the good agreement already observed for vol- 
umes |g, |9| , we tested this functional form to daily data 
of volatilities reported in the literature 0. The results 
are shown in Fig. |3 In the models we are generalizing, 
the variable x is usually identified with the variance or 



squared volatility (x 



Then, the resulting PDF for 



v is 



P(v) 



Po 



{v/v ) 



2i/+l 



[l + (q-l)(v/v o y 



(20) 



with p = 2(2-g)"(«-i)r( 5 i T -i)/Kr( 1 ,+i)r( ? i T - 1 /-i)]. 

Fig. |3 shows an excellent agreement between theoreti- 
cal and empirical PDFs, for the full range of data. Notice 
that the very central part of the distribution is parabolic 
in the log- log plot, then a poor statistics at the tails may 
mislead to think that the distribution is log-normal. 
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FIG. 3: PDF of normalized volatility of stocks traded in US 
equity market (data from Ref. The full line corresponds 
to a fitting by the theoretical PDF given by expression 1201 , 
Fitting parameters are (q, v,v„) ~ (1.178, 2.20, 0.097). Insets: 
linear-linear and log-log representation of the same data. 



Underlying dynamics 

The satisfactory agreement between the empirical data 
and Eq. (|2*U|> suggests that processes similar to either A 
or B may rule squared-volatility evolution. Hence, let 
us look at the explicit form of the ILEs associated to 
processes A and B: 



dx = —j\x — 01 — dt + iiyfxdz + adw . (21) 
x 



B: dx = — ^[x — 6]^rdt + p—=dz 
x \ x 



adw . 



(22) 



The first term in each ILE represents the deterministic 
restoring force with respect to the level 9. It derives 
from a confining potential of the form j[x — 8lnx] (A) 
or 7 [In a; + 9/x] (B). In both cases, the potential has a 
minimum located at x — 9 and is divergent at x = 0. 
Average values are 



Wa = 



g + l/7a 

1 - 1/7m 



and 



Wb 



1 - 2/ 7 o 



(23) 



both averages are greater than 9 and coincide only in the 
limit of relatively small noise amplitudes (7 >> a 2 ,/i 2 ). 
Moments (x n ) are finite only if 7^ > n (A) or j a > n + 1 
(B). In particular, the second moment is 



<* 2 



<* 2 



i a {9 + ih a ){e + 2/ la ) 

7^(1-1/7^(1-2/^) 

7m(A + 2/7/0(0 + 3/7m) 
7o(l - 2/ 7q )(1 - 3/ 7 „) 



(24) 



(25) 



In model A, increasing(decreasing) amplitude of the ad- 
ditive (multiplicative) noise, increases the width of the 
distribution, whereas model B presents opposite behav- 
ior. Thus, for instance, the additive noise has a confining 
effect in process A, opp osite to the effect observed in 
processes with null fl|l8j|. 

On the other hand, the distribution has a maximum 

at 



1 + 1/7, 



and x% ax 



1 

7m 



(26) 



Notice that the additive noise does not affect the mode, 
as expected. The most probable value of distribution A 
shifts to the right with increasing multiplicative ampli- 
tude, while in distribution B the opposite tendency oc- 



curs. From Eqs. and ||2EJ, 



< 9a < (x)a, while 



9b < x r £ ax < (x)b- That is, in model A, the reference 
value 9 represents a typical value comprised between two 
central measures, which does not hold in model B. This 
observation, in addition to the positivity of s, point to 
model A as a more realistic long-term process. 

The fitting parameters in Fig.|31lead to (0,7p,7 q )a — 
(0.50,2.4,4.4) or (9, 7^ 7ct ) B ^ (0.19,6.3,3.4). In both 
cases, 7 m ,7q > 1, as expected for regulated markets. 
While (v) — 1, because empirical volatility is normal- 
ized, (x) = (v 2 ) ~ 1.3 and the mode is x max ~ 0.35, 
consistently with Eqs. P5 |l -l(2l) jl , 

Numerical i nteg ration of ILEs (|21l) and (1221) , by stan- 
dard methods |la|. shows that both processes produce 
time series with bursting or clustering effects, as observed 
in real sequences. However, process B may present, for 
some values of the parameters, a kind of ergodicity break- 
ing, with large jumps to a state basically governed by 
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additive noise. This occurs because, once x jumps to a 
high value, both the restoring force and the effective am- 
plitude of multiplicative noise become small as to pull 
x back to its reference level. Then, relaxation is slowed 
down and the regime of high volatility persists for long 
time stretches. Although a process with s < is not 
expected to be a realistic model for very long time inter- 
vals, it can model, for instance, the transient behavior of 
the market around crashes. In fact, process B yields akin 
crises. After the crash occurs, this drastic event might 
switch the system back to a s > regime. 

VI. FINAL REMARKS 

We have analyzed stochastic models of a quite general 
form, with algebraic restoring force and algebraic multi- 
plicative noise. A further generalization with the inclu- 
sion of an extra source of noise, of standard Wiener type, 
has also been analyzed. These additive-multiplicative 
processes are built on the basis of realistic features: The 
multiplicative noise describes innovations generated by 
endogenous mechanisms that amplify or attenuate a ran- 
dom signal, depending on the internal state of the sys- 
tem. Whereas, the additive noise encodes a direct in- 
fluence of external random fields such as news or sponta- 
neous fluctuations due to speculative trading. One of the 
goals of this work was to study systematically the PDF 
asymptotic solutions of these generalized models. We 
have shown that the inclusion of additive noise gives rise 
to new PDF shapes, with a richer spectrum of cross-over 
behaviors and, in particular, two-fold power-law decays. 
The shapes of the PDFs are governed by the effective 
market rules parametrized by r and s. These parameters 
describe the algebraic nature of the global mean-reverting 
strength of the market and the informational coupling 
among the traders, respectively. On the other hand, 
power-law exponents and coefficients of exponential-like 



functions depend also on the reduced parameters 7 M , j a 
and on 9. This means that one may expect universal 
behavior among markets that share similar rules (same 
r and s) and same rescaled restoring parameters, for a 
properly normalized reference level 9. Summarizing, the 
additive-multiplicative processes given by Eq. (|10|l pro- 
vide a general realistic framework to describe the shape 
of empirical distributions for financial, as well as, for 
physical systems. An illustrative application to empirical 
volatility data was presented in Sect. [V] showing excel- 
lent results. 

The statistical description of a market should include 
its dynamical properties such as the temporal decay of 
correlations. In real time series of volatilities^ and 
volumes 0, power-law decaying correlations have been 
observed. It is worth noting that stochastic processes 
with additive-multiplicative structure (without mean re- 
version) are being currently studied in connection with 
a generalization of standard (Boltzmann-Gibbs) statisti- 
cal mechanics, recently proposed by C. Tsallis [25(. The 
PDFs associated to this new formalism generalize the ex- 
ponential weights, namely, e xp q ( — x) = (1 — [1 — q]x) 1 -i 
[entering as a factor in Eq. I|19l) ]. The time series aris- 
ing from additive-multiplicative processes without mean 
reversion present strong correlations that prevent conver- 
gence to either Gauss or Levy limitsj^ and lead to q- 
Gaussian distributions. This suggests that similar corre- 
lations may persist in mean-reverting processes with the 
additive-multiplicative character. Once Eq. (|10fl leads 
to PDFs in such a good agreement with empirical ones, 
it is worth performing a detailed study and comparison 
of real and artificial time series to test the models with 
respect to the dynamics. Elucidating this point deserves 
a careful separate treatment. 
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